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postshot  recalculation  of  the  event.  The  pretest  calculational  grid  extended 
only  to  a  depth  of  about  300  feet.  Examination  of  possible  travel  paths  showed 
that  reflection  from  deeper  interfaces  could  have  arrived  in  time  to  influence 
the  ground  motions  measured  during  Event  III.  Therefore,  it  was  recommended 
that  the  recalculation  grid  extend  to  depths  of  500  to  600  feet,  i.e.,  into  the 
Precambrian  basement.  Tests  on  samples  taken  just  after  the  event  showed  that 
water  content  changes  due  to  the  wet  weather  conditions  at  shot  time  substan¬ 
tially  affected  the  compressibility  of  the  upper  2  to  3  feet  of  overburden  soil. 
Thus,  it  was  recommended  that  the  overburden  soil  be  idealized  into  two  layers 
rather  than  into  one  as  or  ginally  recommended  for  preshot  calculations.  An 
extensive  reanalysis  of  the  available  surface  retraction  survey  data,  along  with 
analysis  of  the  Event  III  ground  motion  data,  resulted  in  revisions  to  the 
entire  seismic  velocity  profile  and,  hence,  to  the  associated  values  of  the 
initial  uniaxial  strain  modulus.  The  stress  levels  associated  with  these  moduli 
were  originally  assumed  to  be  quite  low,  e.g.,  less  than  10  psi  for  the  Kayenta 
Formation  materials.  A  lower  bound  value  of  50  psi  was  recommended  as  being 
more  reasonable  for  the  initial  Kayenta  "precursor’'  stress.  This  phenomenon 
appears  to  be  a  function  of  loading  rate.  Tests  on  the  Kayenta  materials  also 
revealed  that  horisontal-to-vertical  anisotropy  undoubtedly  affected  the  ground 
motion  results,  but  the  problem  of  how  to  utilize  the  various  horizontal  and 
vertical  data  to  specify  meaningful  "effective"  or  "average"  isotropic  property 
values  still  remains.  As  a  result  of  recent  test  data  from  Lawrence  Livermore 
Laboratory,  however,  substantial  changes  have  been  made  in  the  recommended 
failure  envelopes,  which  now  reflect  a  highly  nonlinear  behavior,  including 
significant  strength  increases  at  high  pressure. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  bi  converted  to 
metric  units  as  follows: 


Multiply 


By 


To  Obtain 


feet 

tons 

feet  per  second 

pounds  (mass)  per  cubic 
foot 

pounds  (force)  per  square 
inch 

kips  (force)  per  square 
inch 


0.3048 

907.1846 

0.3048 

16.0185 

0.6894757 

0.6894757 


meters 

kilograms 

meters  per  second 

kilograms  per  cubic  meter 

neurons  per  square  centimeter 

kilonewtons  per  square 
centimeter 


pounds  (force)  per  square 
inch  per  second 


0.6894757  newtons  per  square  centimeter 
per  second 


CHAPTER  1 


INTRODUCTION 


1.1  BACKGROUND 

A  geologic  profile  for  the  Mixed  Company  site  and  associated  con¬ 
stitutive  properties  for  each  layer  in  this  profile  were  developed  dur¬ 
ing  FY  1972  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES).  The  idealized  profile  and  properties  (Reference  l)  were  dissem¬ 
inated  in  June  1972  for  use  in  developing  fits  to  constitutive  models 
that  were  in  turn  used  for  pretest  ground  motion  calculations  (Refer¬ 
ence  2)  of  Event  III,  a  500-ton*'  high  explosive  test.  Since  these  prop¬ 
erties  were  disseminated,  additional  field  and  laboratory  test  data,  as 
well  as  ground  shock  data  from  the  event  itself  and  from  the  Air  Force 
Weapons  Laboratory's  (AFWL)  Cylindrical  In  Situ  Test  (CIST)  explosive 
experiment,  have  become  available. 

Ground  shock  data  obtained  both  by  WES  and  AFWL  were  presented  at 
the  Mixed  Company  project  review  meeting  in  March  1973.  The  pretest 
ground  motion  predictions  did  not  agree  very  well  with  either  set  of 
measurements.  The  calculated  arrival  times  for  the  initial  ground  mo¬ 
tions  were  much  later  than  those  indicated  by  the  field  measurements , 
and  the  general  characteristics  of  the  calculated  wave  forms  did  not 
match  those  of  the  field  data.  A  number  of  assumptions  were  made  in  the 
process  of  developing  the  preshot  profile  and  property  idealizations 
(Reference  3),  which  if  not  valid  could  well  account  for  the  observed 
discrepancies  between  the  field  measurements  and  the  calculated  ground 
motions.  It  should  also  be  noted  that  the  mathematical  models  had  to  be 
hurriedly  fit  in  order  to  meet  the  schedule  for  the  pretest  calculation, 
and  they  did  not  closely  replicate  all  facets  of  the  recommended  consti¬ 
tutive  properties.  Also,  the  calculational  grid  did  not  extend  deep 


A  table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  9. 
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enough  to  incorporate  the  reflecting  interface  between  the  Triassic  sed¬ 
imentary  materials  and  the  Precamibrian  basement  rock.  Thus,  there  are 
sufficient  uncertainties  to  warrant:  (l)  a  reevaluation  of  the  profile 
and  properties  in  the  light  of  the  newly  available  data,  and  (2)  a  re¬ 
calculation  of  Event  III  with  a  much  morfe  precise  fit  to  the  constitu¬ 
tive  properties  and  profile  resulting  from  this  reevaluation. 

1.2  PURPOSES 

The  purposes  of  this  study  were  to  evaluate  the  validity  of  the 
preshot  Mixed  Company  Event  III  site  profile  and  the  associated  set  of 
constitutive  properties  in  light  of  new  data  and  to  make  appropriate 
recommendations  for  a  revised  site  profile  and  a  revised  set  of  consti¬ 
tutive  properties  for  use  in  a  postshot  recalculation  of  the  event. 

1.3  SCOPE 

This  paper  presents  the  results  of  a  reevaluation  of  the  Mixed  Com- 

'  '  '/ 

pany  site  profile  and  constitutive  properties  utfili-zing  both  field  and 
laboratory  data.  Revisions  to  the  subsurface  profile  were  based  pri¬ 
marily  on  the  water  contents  determined  in  the  soil  overburden  at  shot 
time.  An  extensive  reanalysis  of  the  seismic  velocity  data  was  con¬ 
ducted.  New  constitutive  property  recommendations  for  the  overburden 
soil  were  developed  based  on  laboratory  tests  conducted  on  specimens 
taken  in  the  field  just  a  few  days  after  the  event.  The  properties  of 
the  Kayenta  Formation  were  revised  to  consider  the  effect  of  loading 
rate  on  the  initial  portion  of  the  uniaxial  strain  (UX)  test  stress- 
strain  curves,  the  effects  of  anisotropy  on  the  observed  triaxial  shear 
(TX)  and  UX  test  behaviors,  and  the  correlation  observed  between  shear 
strength  obtained  from  laboratory  tests  and  strength  obtained  from  an 
in  situ  test. 
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CHAPTER  2 


SITE  PROFILE 

In  developing  the  pretest  calculation  profile  for  the  Mixed  Company 
Event  III  site,  many  assumptions  „e.re  necessary  in  order  to  translate 
the  available  data  into  an  idealized,  horizontally  layered  profile.  The 
purpose  .f  this  chapter  is  to  evaluate  these  assumptions  in  the  light  of 
additional  information  and  analyses  and  to  recommend  appropriate  revi¬ 
sions.  As  wa '»  the  cate  for  the  original  profile,  the  revised  profile  is 
intended  to  be  an  idealized  representation  of  conditions  along  the  WES 
gage  line.1 

2.1  DEPTHS  TO  MAJOR  GEOLOGIC  INTERFACES 

Four  major  geologic  interfaces  were  originally  defined  for  the  site 
profile.  The  first,  at  a  depth  of  5  feet,  depicted  the  change  from  soil 
overburden  to  the  complex  Kayenta  Formation  of  intermixed  siltstones, 
sandstones,  mudstones,  and  conglomerates.  At  70  feet,  the  sandstones  of 
the  Wingate  Formation  were  assumed  to  be  encountered,  followed  by  Chinle 
Formation  siltstone  at  400  feet.  Finally,  an  interface  with  the  Preeam- 
brian  basement  complex  was  set  at  a  depth  of  500  feet.  Nc  reason  has 
been  found  to  alter  these  interface  depths.  The  pretest  calculational 
grid,  however,  extended  only  to  a  depth  of  approximately  300  feet  and 
thus  included  only  the  overburden  soil,  the  Kayenta  Formation,  and  a 
portion  of  the  Wingate .  Since  that  time ,  several  calculations  have  been 
made  for  other  piojects  that  indicate  that  near-surface  ground  motions 
can  be  significantly  affected  by  reflections  from  very  deep  geologic  in¬ 
terfaces.  This  determination  prompted  an  examination  of  the  profile  to 
see  if  signals  from  the  omitted  Mixed  Company  interfaces  were  being 


Two  gage  lines  were  installed  for  Event  III,  one  by  VJES  and  one  by 
AFWL.  A  rolled  fill  was  constructed  along  the  WES  line  in  an  attempt 
bo  provide  a  uniform  5--foot-thick  layer  of  soil  overburden;  fill  was 
not  placed  along  the  AFWL  line,  and  the  natural  soil  cover  varied  in 
thickness  from  3-1/2  to  5  feet. 
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recorded.  Figure  1  shows  the  subsurface  profile  to  a  range  of 
1,^400  feet  end  to  a  depth  of  300  feet.  Also  shown  are  travel  paths  that 
are  possible  for  some  of  the  reflected  waves  and  the  surface  arrival 
times  of  these  waves  based  on  an  assumed  7*500  ft/sec  wave  speed  in  both 
the  Kayenta  and  Wingate  Formations.  This  approximate  model  shows  that 
signals  reflected  from  the  Chinle  and  Precambrian  layers  could  have  ar¬ 
rived  in  time  to  influence  the  ground  motions  measured  during  Event  III, 
since  significant  near-surface  motion  occurred  for  at  least  500  msec  at 
most  of  the  ranges  shown  in  the  figure.  It  is  therefore  strongly  recom¬ 
mended  that  the  grid  used  in  any  Mixed  Company  recalculation  extend  into 
the  Precambrian  basement  complex. 

2.2  IDEALIZATION  OF  PRIMARY  GEOLOGIC  UNITS 

As  denoted  above,  the  site  profile  encompasses  five  primary  geo¬ 
logic  units,  i.e.,  the  soil  overburden,  Kayenta  sandstones,  Wingate 
sandstones,  Chinle  siltston ;s ,  and  the  Precambrian  basement  rocks.  In 
the  idealized  pretest  profile,  the  overburden  soil  was  represented  as  a 
single  horizontal  layer,  Layer  I.  The  underlying  Kayenta  Formation, 
however,  was  subdivided  into  four  layers:  Layer  II,  representing  the 
weathered  siltstcne;  Layer  III,  the  upper  one-fourth  of  the  unweathered 
Kayenta  rock;  Layer  IV,  an  artificial  layer  representing  the  softer 
Kayenta  material  occurring  ’* mdomly  throughout  the  formation;  and 
Layer  V,  the  lower  three-fourths  of  the  unweathered  Kayenta  rock.  The 
Wingate,  Chinle,  and  Precambrian  Formations  were  each  represented  by  one 
layer,  Layers  VI,  VII,  and  VIII,  respectively.  In  making  this  idealiza¬ 
tion,  it  was  assumed  that  the  water  content  and  density  of  the  con¬ 
structed  berm  would  match  those  of  the  natural  soil  ov^fburden  and  also 
that  the  water  content  and  other  compositional  properties  found  to  exist 
at  the  time  samples  were  extracted  from  the  site  would  be  the  same  as 
those  existing  at  the  time  of  actual  detonation. 

Wet  weather  conditions  existed  at  the  site  just  prior  to  Event  III, 
and  the  water  content  of  the  upper  2  to  3  feet  of  overburden  soil  in¬ 
creased  from  an  average  of  about  7  percent,  as  determined  at  the  time 
preshot  properties  were  developed,  to  about  15  percent  at  the  time  of 
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the  shot.  The  variations  in  water  content  and  dry  density  with  depth  in 
the  top  5  feet  in  a  section  in  natural  soil  are  shown  in  Figure  2. 

From  the  surface  to  a  depth  of  approximately  2  feet,  the  water  content 
at  shot  time  appears  to  have  been  approximately  15  to  1 6  percent;  at 
depths  below  3  feet,  however,  the  available  data  indicate  that  the  water 
content  may  not  have  changed  significantly  from  the  original  value  of 
about  7  percent.  The  variation  in  dry  density  with  depth  shown  in  Fig¬ 
ure  2  indicates  that  no  significant  change  in  that  composition  prop¬ 
erty  had  occurred  since  the  time  of  the  pretest  investigation.  Fig¬ 
ure  3  shows  the  dry  density  and  water  content  profile  for  a  composite 
section  of  fill  and  natural  soil.  The  water  content  profile  is  essen¬ 
tially  the  same  as  that  shown  in  Figure  2;  the  dry  density  profile, 
however,  is  quite  different.  The  density  of  the  fill  itself  appears  to 
be  very  uniform,  but  the  upper  1/2  to  1  foot  of  natural  soil  underlying 
the  fill  apparently  increased  in  density  due  to  compaction  that  occurred 
during  construction  of  the  berm. 

Primarily  on  the  basis  of  the  data  shown  in  Figures  2  and  3,  it 
is  recommended  that  the  soil  overburden  along  the  main  WES  gage  line  be 
idealized  into  two  horizontal  layers  for  postshot  calculations  rather 
than  into  one  as  originally  recommended  for  preshot  calculations. 

The  postshot  water  content  and  densities  of  the  Kayenta  materials 
did  not  differ  from  the  pretest  values  sufficiently  to  warrant  a  revi¬ 
sion  in  the  subdivision  of  the  Kayenta  Formation.  Although  the  wet  den¬ 
sity  of  the  upper  one-fourth  of  the  unweathered  Kayenta  material  was 
lower  by  approximately  2  pcf  than  the  pretest  density ,  it  was  assumed 
that  the  new  value  simply  represented  a  better  definition  of  the  density 
for  that  material  rather  than  an  alteration  due  to  the  blast  or  change 
in  climate.  None  of  the  clayey  conglomerate  material,  which  was  arbi¬ 
trarily  placed  at  the  18-  to  22-foot  depth  to  form  Layer  IV  in  the  pre¬ 
test  idealized  profile,  was  encountered  in  any  of  the  four  20-ft-deep 
posttest  borings.  Therefore,  the  composition  properties  for  this  mate¬ 
rial  could  not  be  reevaluated. 

No  additional  new  data  were  obtained  in  order  to  evaluate  the  pre¬ 
shot  idealizations  recommended  for  the  deep  formations  underlying  the 
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Kayenta.  However,  the  geologic  information,  the  limited  dry  density 
data  obtained  on  field  samples  taken  at  outcrops,  and  the  water  content 
values  that  were  previously  ass’imed  were  reviewed.  As  a  result,  it  was 
again  decided  that  the  Wingate,  Chinle,  and  Precambrian  Formations 
should  not  be  subdivided  but  that  each  can  be  adequately  depicted  as  a 
single  idealized  horizontal  layer  for  caleulational  purposes. 

2.3  RECOMMENDED  POSTSHOT  PROFILE 

The  recommended  postshot  profile  and  composition  properties  are 
summarized  in  Table  1,  which  lists  the  numb^is  of  the  old  preshot 
layers  and  those  of  the  new  postshot  layers,  the  material  description, 
the  depth  range  for  each  layer,  and  the  values  of  wet  density,  water 
content,  and  the  volume  of  air  now  judged  to  most  nearly  represent  con¬ 
ditions  along  the  WES  gage  line.  There  are  two  significant  differences 
between  this  profile  and  the  pretest  profile.  First,  the  soil  overbur¬ 
den  has  been  subdivided  into  two  layers,  Layers  1A  and  IB.  However,  if 
the  caleulational  grid  cannot  be  zoned  fine  enough  to  effect  this  recom¬ 
mended  subdivision,  then  a  composite  layer.  Layer  1C,  has  also  been  de¬ 
fined  to  represent  the  total  zone  of  fill  and  natural  soil  overburden. 

Second,  it  is  recommended  that  the  artificial  Layer  IV,  formerly 
located  between  a  depth  of  18  and  22  feet,  be  eliminated  from  the  post¬ 
shot  calculation  profile,  at  least  initially.  With  the  uncertainties 
surrounding  the  definition  of  seismic  velocities  and  constitutive  prop¬ 
erties,  as  will  be  discussed  later,  it  appears  that  the  most  practical 
approach  would  be  to  keep  the  profile  as  simple  as  possible  for  the 
planned  recalculation.  However,  the  basic  assumption  of  lumping  small 
zones  of  material,  which  cannot  be  incorj  crated  in  the  profile  due  to 
computational  grid-size  limitations,  into  one  or  more  large  layers 
should  be  evaluated  in  future  parametric  studies. 
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Figure  1  Subsurface  profile  showing  possible  paths  of  waves 
reflected  from  deep  geologic  interfaces  and  corresponding 
arrival  times  at  the  ground  surface. 
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Figure  Profile  of  natural  soil  overburden  layer 
shoving  variation  of  water  content  and  dry  density 
with  depth. 
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Figure  3  I  roiile  of  composite  '.'ill  and  natural  soil 
overburden  layer  showing  variation  of  water  content 
and  dry  density  with  depth. 
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CHAPTER  3 


SEISMIC  VELOCITIES 

If  the  seismic  velocities  used  with  the  calculation  profile  do  not 
match  those  of  the  real  site  profile,  then  calculated  arrival  times  ob¬ 
viously  cannot  agree  with  field-measured  data.  Neither  can  the  calcu¬ 
lated  wave  forms  match  the  field  wave  forms,  since  the  calculated  wave 
arrivals  will  be  incorrectly  phased.  For  this  reason,  an  important  key 
to  an  adequate  calculation  of  Mixed  Company  Event  III  is  the  correct 
definition  of  seismic  velocities.  The  analyses  originally  applied  to 
the  seismic  velocity  data  from  22  surface  refraction  lines  are  docu¬ 
mented  in  the  summary  report  (Reference  1);  however,  these  data,  along 
vrith  those  available  from  four  other  refraction  lines  surveyed  near  the 
test  site,  have  now  been  reexamined  with  respect  to  the  overall  topog¬ 
raphy  of  the  area.  In  addition,  da '.a  from  the  AFWL  CIST  experiment  have 
provided  new  information  on  horizontal  wave  velocities,  and  the  arrival¬ 
time  data  from  Event  III  itself  have  been  analyzed  as  a  refraction  sur¬ 
vey  to  deduce  possible  seismic  velocity  profiles.  The  purpose  of  this 
chapter  is  to  present  the  results  of  these  various  examinations  and  to 
recommend  seismic  velocity  values  for  the  site  profile  developed  in 
Chapter  2  for  postshot  calculations. 

3.1  REANALYSIS  OF  SURFACE  REFRACTION  SURVEY  DATA 

A  plan  view  of  the  Mixed  Company  site  showing  seismic  survey  loca¬ 
tions  with  respect  to  the  Event  III  ground  zero  (GZ),  the  WES  and  AFWL 
gage  lines,  and  the  AFWL  CIST  experiment  is  shown  in  Figure  H.  Seis¬ 
mic  velocity  data  from  21  hammer-source  surface  refraction  lines,  each 
150  feet  long  with  15-foot  geophone  spacings ,  and  one  570-foot-long  ex¬ 
plosive  line  wc^e  obtained  from  the  immediate  vicinity  of  GZ  (Refer¬ 
ence  b) .  The  short-line  (150-foot-long)  data  indicated  velocities  rang¬ 
ing  from  2,100  to  3,200  ft/sec  in  the  upper  9  feet  (i.e.,  a  combination 
of  soil  overburden  and  weathered  Kaye  its  siltstone )  and  velocities  rang¬ 
ing  from  6,000  to  8,800  ft/sec  in  bhe  immediately  underlying  unweathered 
Kayenta  rock.  In  the  pretest  analysis,  a  value  of  1,800  ■f'c/sec  was 


19 


estimated  for  the  overburden  soil,  a  value  of  2,500  ffc/sec  for  the 
weathered  Kayenta,  and  values  of  7,100  ft/sec  for  the  upper  one-fourth 
of  the  unweathered  Kayenta  and  7? 500  ft/sec  for  the  lower  three-fourths. 

The  long-line  (570-foot-long)  survey  indicated  a  zone  of 
lH ,000-ft/ sec  material  at  a  depth  of  about '100  feet,  but  this  value  was 
not  recommended  for  preshot  calculations  since  it  was  believed  to  be  too 
high  to  be  representative  of  the  low-density  Wingate  sandstone  expected 
at  that  depth;  a  pretest  value  of  7,4-00  ft/sec  was  recommended.  None 
of  the  refraction  lines  were  long  enough  to  provide  seismic  Vv "  ocit-ies 
for  the  Chinle  siltstone  and  the  Erecambrian  basement;  values  of  9,000 
and  15,000  ft/sec,  respectively,  were  estimated  for  these  final  two 
layers  of  the  pretest  profile. 

Seismic  surveys,  each  consisting  of  25-,  50- ,  and  600-foot~long 
refractions  lines ,  were  also  run  at  the  four  locations  shown  in  Figure  4 
that  lie  between  the  main  WES  gage  line  and  the  northwest-trending  mesa. 
The  data  from  these  surveys  were  not  considered  in  the  initial  analysis; 
however,  they  were  subsequently  examined  since  the  short  lines  -»  .re  run 
with  5-foot  geophone  spacings  that  enabled  the  seismic  velocities  in  the 
overburden  soil  (ranging  from  1,00.0  to  1,400  ft/sec  with  an  average  of 
1,100  ft/sec)  to  be  distinguished  from  those  in  the  weathered  siltstone 
(ranging  from  1,800  to  5,000  ft/sec  with  an  average  of  3,300  ft/sec). 

Velocities  in  the  deeper,  more  competent  Kayenta  materia.,,  seemed 
to  be  related  to  the  areal  topography.  Four  cross  sections  through  the 
site  are  depicted  in  Figure  5*  The  geologic  formation  interfaces  were 
estimated  from  visual  observations  of  the  various  mesa  and  canyon  walls; 
the  Kayenta-Wingate  interface  appeared  to  roughly  follow  the  surface 
slope.  The  available  seismic  survey  data  have  >,een  projected  onto  two 
of  the  cross  sections  (Sections  C-C  and  D-D);  the  lower  velocity  zones 
(4,000  to  6,000  ft/sec)  of  the  upper  Kayenta  material  appeared  to  be 
confined  lx>  the  lower  elevation  drainage  patterns,  while  higher  veloc¬ 
ities  (7, C00  to  8,000  ft/sec)  indicative  of  a  more  resistant  material 
were  recorded  in  the  higher  elevation  area  around  GZ.  Velocities  for 
the  lower  Kayenta  materials  in  either  case  appeared  to  be  on  the  order 
of  8,000  to  10,000  ft/sec.  As  shown  in  Figure  4,  the  WES  gage  lines 
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are  located  entirely  along  a  knoll  where  75000-  to  8,000-ft/sec  veloc¬ 
ities  should  typify  the  upper  unweathered  Kayenta  material.  Along  the 
AFWL  gage  line,  however,  two  profiles  may  exisc.  Close  to  GZ 
(<300  feet),  the  profile  should  "be  similar  to  that  under  the  WES  gage 
line,  but  farther  out  as  the  terraii.  drops  off  the  upper  Kayenta  veloc¬ 
ities  may  be  on  the  order  of  4,000  to  6,000  ft/sec. 

Although  the  four  early  600-foot  line  surveys  easily  penetrated  the 
Wingate  formation,  none  detected  a  velocity  increase  over  the  lower 
Kayenta.  The  Wingate  could  have  had  a  lower  seismic  veiocity  than  the 
Kayenta,  but  this  difference  could  only  have  been  detected  with  an  up- 
hole  rather  than  a  refraction- type  field  survey.  Seismic  velocity  mea¬ 
surements  made  in  the  laboratory  on  samples  of  Wingate  sandstone  were  in 
fact  lower  than  similar  measurements  made  on  the  overlying  Kayenta  sam¬ 
ples.  The  four  lines  were  too  short,  however,  to  penetrate  the  Chinle 
and  Precambrian  Formations;  however,  laboratory  tests  have  indicated 
that  seismic  velocities  of  approximately  11,000  and  18,000  ft/sec,  re¬ 
spectively,  should  be  expected. 

3.2  ANALYSIS  OF  EXPLOSIVE  TEST  DATA 

On  28  September  1972,  AFWL  conducted  a  35-foot-deep  CIST  experiment 
about  1  mile  to  the  southeast  of  the  Mixed  Company  Event  III  GZ  (see 
Figure  4.).  The  first  arrival-time  data  indicated  horizontal  seismic 
velocities  in  the  upper  unweathered  Kayenta  material,  ranging  from  ab^ut 
10,000  ft/sec  to  about  20,000  ft/sec  (Reference  5).  These  velocities 
were  much  higher  than  those  indicated  by  the  surface  refraction  surveys. 
Since  they  were  measured  along  a  horizontal  plane,  they  may  indicate  a 
significant  horizontal-to-vertical  anisotropy  for  the  Kayenta  material. 
Even  so,  20,000-ft/sec  velocities  are  generally  associated  with  sound 
igneous  rocks,  such  as  granite,  rather  than  with  relatively  soft  sedimen¬ 
tary  rocks  such  as  siltstone  and  sandstone.  A  cross-hole  seismic  survey 
at  the  Mixed  Company  site  is  definitely  recommended  to  assist  in  resolv¬ 
ing  this  uncertainty.  Since  the  CIST  measurements  were  obtained  at 
higher  elevations  than  any  of  the  other  measurements ,  it  may  also  be 
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possible  that  the  velocity  differences  are  due  to  material  and/or  geo¬ 
logic  differences. 

Event  III  was  detonated  on  13  November  1972.  Shock-front  diagrams 
were  constructed  by  both  WES  (Reference  6)  and  .\FWL  (Reference  7)  based 
on  the  arrival  times  recorded  by  their  respective  motion  gages.  These 
diagrams  indicated  that  the  soil  overburden  had  an  average  seismic  veloc¬ 
ity  of  about  1,3^0  ft/sec;  this  agrees  quite  well  with  the  1,100-ft/sec 
velocity  value  determined  by  the  refraction  surveys,  ".'he  date  indicate 
velocities  on  the  order  of  3,500  ft/sec  for  the  weathered  siltstone 
stratum;  again,  this  is  in  excellent  agreement  with  the  3,300-ft/sec 
average  determined  from  the  refraction  data.  The  data  also  indicate 
velocities  for  the  unweathered  Kay^nta  materials  in  the  8,000-  to 
10,000-ft/see  range,  which  are  also  consistent  with  the  refraction 
survey  results. 

The  only  gages  placed  deep  enough  to  be  in  the  Wingate  sandstone 
were  those  directly  under  GZ.  It  is  possible  that,  in  this  region  where 
the  travel  paths  es.-entially  paralleled  the  gage  columns,  the  grout  used 
to  fill  uhe  gage  columns  may  have  influenced  the  arrival  times,  since  it 
had  a  higher  P-w.ave  velocity  (i.e.,  about  10,000  ft/sec)  than  the  over¬ 
all  mass  of  natural  materials.  However,  the  shock-front  diagram  data 
are  insufficient  to  permit  a  definite  conclusion  in  this  regard. 

The  near-surface  (l-l/2-foot-deep)  motion  gage  results,  however, 
can  be  analyzed  as  a  long-line,  surface  refraction  survey  to  extract 
some  seismic  velocity  information  for  the  Wingate  and  deeper  forma¬ 
tions.1  One  such  analysis  of  the  data  obtained  along- the  WES  gage  lines 
(LN302  and  LN311W)  was  reported  by  Ballard  and  Leach  (Reference  8).  It 
indicated  a  single  570-foot-thick  stratum  of  8,300-ft/sec  material  over- 
lying  l8,000-ft/sec  material  (see  Example  A  in  Figure  6).  Other  in¬ 
terpretations  or  fits  to  the  data  are  possible,  as  indicated  by 


This  technique  will  not  provide  valid  information  for  the  materials 
above  a  depth  of  about  60  to  70  feet  since  the  airblast-dominated 
range  of  superseismic  motion  is  the  "effective"  geophone  spacing  for 
the  survey. 
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Examples  B,  C,  and  D  in  the  insert  in  Figure  ( .  When  the  layer  veloc¬ 
ities  do  not  monotonic ally  increase  with  depth,  as  is  bhe  case  with  Ex¬ 
amples  B  and  C,  the  layer  interface  depths  cannot  he  calculated  from  the 
refraction  equations;  therefore,  only  a  i^agh  estimate  is  possible. 

3.3  RECOMMENDED  POSTSHOT  SEISMIC  VELOCITIES 

The  shock-front  or  arrival -time  diagrams  from  Event  III  and  the  re¬ 
cently  analyzed  short-line  refraction  surveys  both  indicate  that  the 
1,800-ft/sec  seismic  velocity  originally  recommended  for  the  soil  over¬ 
burden  was  too  high;  a  value  of  1,300  ft/sec  Is  now  recommended  for 
postshot  calculation  models.  Conversely,  the  arrival-time  data  and  the 
short-line  refraction  data  indicated  that  the  2,500-ft/sec  velocity  pre¬ 
viously  recommended  for  the  immediate?^y  underlying  weathered  silts  tone 
was  too  low;  a  value  of  3,500  ft/sec  is  now  recommended. 

Reanalysis  of  all  of  the  available  surface  refraction  data  indi¬ 
cates  that  velocities  in  the  unweathered  Kayenta  materials  depend  on  the 
areal  topography,  with  values  within  the  ranges  shown  in  Figure  7* 

Along  the  WES  gage  line,  a  specific  value  of  8,000  ft/sec  is  now  recom¬ 
mended  for  the  upper  one-fourth  of  the  unweathered  Kayenta,  and 
10,000  ft/sec  is  recommended  for  the  lower  three-fourths.  Previously 
recommended  values  were  7,1°°  and  7,500  ft/sec,  respectively. 

As  shown  in  Figure  7,  a  variety  of  field  values  were  obtained 
for  depths  associated  with  the  Wingate  sandstone,  but  all  exceeded 
8,000  ft/sec;  laboratory  tests,  however,  gave  values  consistently  lower 
than  8,000  ft/sec.  A  velocity  of  8,000  ft /sec  is  now  recommended  in 
lieu  of  the  previous  value  of  7, **00  ft/sec.  Based  solely  on  laboratory 
test  data,  Chinle  siltstone  velocities  are  now  believed  to  be  on  the 
order  of  11,000  rather  than  9,000  ft/sec  as  originally  recommended. 
Finally,  the  Precambrian  basement  velocity  is  now  believed  to  be  at 
least  18,000  ft/sec  rather  than  15,000  ft/sec.  This  conclusion  was 
derived  from  both  the  Event  III  arrival -time  analysis  and  the  recent 
laboratory  test  results. 

These  current  best  estimates  of  the  seismic  velocities  associated 
with  each  of  the  layers  in  the  idealized  geologic  profile  beneath  the 
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WES  gage  line  are  summarized  in  Table  2.  Previously  recommended  values 
are  also  tubulated  for  comparison. 
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Figure  4  Plan  view  of  Mixed  Company  site  showing  seismic  survey 
locations  with  respect  to  the  Event  III  GZ,  the  WES  and  AFWL  gage 
lines,  and  the  juWL  CIST  experiment. 
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Figure  5  Cross  sections  through  Mixed  Company  site  (see 
Figure  *i)  indicating  seismic  velocities  and  estimated 
seismic  and  geologic  interfaces. 
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NOTE:  NUMBERS  BY  CURVES  DENOTE 
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Figure  7  Summary  plot  of  seismic  velocity  profiles  in  Mixed  Company  rocks  (based  on  interpretations 
of  various  data)  compared  with  estimated  geologic  profile. 


CHAPTER  it 


CONSTITUTIVE  PROPERTIES 

The  constitutive  or  mechanical  properties  recommended  for  use  in 
preshot  ground  motion  predictions  were  selected  so  as  to  oe  compatible 
with  other  assumed  properties,  such  as  water  content  and  seismic  veloc¬ 
ity.  Obviously  now,  if  the  original  assumptions  regarding  water  con¬ 
tents  and  seismic  velocities  have  been  modified,  the  associated 
stress-strain  and  shear  strength  relationships  must  be  modified.  The 
mathematical  constitutive  model  formulations  employed  in  the  preshot 
calculations  were  based  on  the  assumption  of  loading-rate-independent 
and  isotropic  material  behavior.  It  was  also  assumed  that  the  labora¬ 
tory  test  results  were  representative  of  in  situ  properties,  and,  in 
some  cases ,  extrapolations  of  the  available  data  were  necessary  in  order 
to  specify  behavior  at  pressures  above  the  range  of  laboratory  test 
values.  The  purposes  of  this  chapter,  therefore,  are  to  examine  new 
data  that  bear  on  these  assumptions  and  suggest  possible  modifications 
to  the  preshot  property  recommendations . 

4.1  EFFECT  OF  WATER  CONTENT  CHANGES 

At  the  time  of  Event  III,  the  surface  soil  at  the  site  was  influ¬ 
enced  by  wet  weather  and  construction  conditions.  As  a  result,  it  was 
recommended  in  Chapter  2  that  the  5-foot-thick  soil  overburden  layer 
originally  specified  be  subdivided  into  two  layers.  The  water  content 
of  the  upper  soil  layer  was  estimated  to  be  approximately  15  percent, 
while  that  of  the  lower  layer  was  estimated  to  be  about  7  percent.  A 
series  of  UX  tests  were  conducted  on  both  undisturbed  and  remolded  spec¬ 
imens  obtained  from  the  site  on  the  day  after  the  event.  As  expected 
(Reference  9),  those  specimens  with  15  percent  water  content  were  ini¬ 
tially  more  compressible  than  those  with  7  p- ’'cent  water  content,  i.e., 
their  initial  UX  response  was  probably  dominates  intergranular  fric¬ 
tion,  which  increases  as  water  content  decreases.  On  the  other  hand, 
the  specimens  with  increased  water  content  contained  fewer  air  voids 
than  did  the  drier  specimens,  which  caused  them  to  stiffen  or  "lock  up" 
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at  smaller  strains  as  these  voids  closed  under  the  applied  stresses. 
Figure  8  shows  a  plot  of  the  recommended  UX  relations  developed  for 
the  two  soil  layers  (Layers  1A  and  IB) ;  an  alternate  layer  (Layer  1C) 
for  the  composite  0-  to  5-foot  depth  is  also  given  for  use  in  the  event 
that  the  calculational  grid  cannot  be  designed  to  accept  the  two-layer 
representation. 

The  TX  failure  envelopes  were  significantly  affected  by  the  changes 
in  water  content.  Test  data  from  remolded  specimens  with  15  percent 
water  content  gave  failure  envelopes  that  achieved  limiting  or  maximum 
principal  stress  difference  values  on  the  order  of  300  to  500  psi  under 
mean  normal  stresses  greater  than  about  500  psi.  The  failure  envelopes 
derived  from  tests  on  specimens  with  only  7  percent  water  content  indi¬ 
cated  a  6,500-psi  limiting  value  for  applied  pressures  on  the  order  of  a 
kilobar  or  greater.  New  yield  surfaces  and  UX  stress  paths  were  devel¬ 
oped  for  the  two  soil  layers  as  a  result  of  these  changes  in  the  shear 
strength  profile. 

k.2  EFFECT  OF  SEISMIC  VELOCITY  CHANGES 

The  initial  uniaxial  strain  modulus  M.  is  assumed  tc  be  directly 

related  to  the  seismic  P-wave  velocity  vp  and  the  mass  density  p  in 

2  V 

the  relation  M.  =  pv_  .  Values  of  M.  were  thus  recalculated  for  all 
l  P  x 

of  the  site  materials  using  the  revised  seismic  velocity  and  density 
values  given  in  Tables  2  and  1,  respectively.  Because  all  of  the  dy¬ 
namic  UX  test  data  indicated  relatively  low  or  soft  initial  moduli,  it 
was  assumed  that  the  extremely  stiff  seismic  moduli  applied  only  for 
very  low  stress  levels  that  were  below  the  resolution  of  the  UX  test 
measurements,  e.g,,  less  than  10  psi  for  the  Kayenta  Formation  materials. 
Although  analyses  of  the  data  are  still  incomplete,  the  general  consen¬ 
sus  at  the  Mixed  Company  Project  Review  Meeting,  held  in  March  1973,  was 
that  the  maximum  amplitudes  of  signals  traveling  with  seismic  velocity 
through  this  formation  in  both  the  AFWL  CIST  experiment  and  Event  III 
were  significantly  larger,  i.e.,  on  the  order  of  200  psi. 

Whereas  the  data  from  dynamic  UX  tests  on  unweathered  Kayenta 
sandstones  did  not  reveal  a  significant  region  of  stiff,  seismic 
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velocity-associated  moduli,  data  from  high-pressure  ( 3-kbar  and  above) 
gun  tests  conducted  by  both  WES  and  the  Stanford  Research  Institute  (SRI) 
certainly  appear  to  indicate  that  moduli  of  this  magnitude  extend  to 
stress  levels  very  much  above  200  psi.  A  comparison  of  the  data  from 
the  two  types  of  tests  is  shown  in  Figure  9*  Assuming  that  both  of 
these  sets  of  laboratory  data  as  well  as  both  sets  of  field  data  are  cor¬ 
rect,  one  possible  explanation  for  the  wide  variation  in  the  stress 
levels  associated  with  th  stiff  initial  loading  moduli  is  that  this 
seismic  "precursor"  variance  is  due  to  loading-rate  effects.  Figure  10 
shows  a  comparison  of  the  loading  times  versus  peak  stress  for  labora¬ 
tory  static  UX  tests,  dynamic  UX  tests,  and  high-pressure  gun  tests  with 
those  estimated  from  the  CIST  and  Event  III  field  data.  The  order  of 
the  data  is  certainly  consistent  with  the  loading-rate  hypothesis. 

To  examine  the  effect  of  loading  rate  on  initial  uniaxial  strain 
response,  plots  of  stress/strain  and  stress/time  from  five  dynamic  UX 
tests  on  similar  unweathered  Kayenta  material  were  enlarged  as  shown  in 
Figure  11.  Three  of  the  five  specimens  had  a  stiffer  initial  response 
than  the  other  two,  but  at  axial  stresses  above  about  20  psi,  the  con¬ 
strained  moduli  from  all  five  specimens  were  essentially  the  same  (Fig¬ 
ure  11a).  Figure  lib  shows  that  +l a  stress  pulses  applied  to  these 
specimens  had  different  initial  loading  rates,  i.e.,  2,000  to 
3,000  psi/sec  for  Pulse  1  as  opposed  to  6,000  to  7,000  psi/sec  for 
Pulse  2\  at  20  psi,  the  loading  rates  for  both  pulse  typos  increased  to 
about  30,000  psi/sec.  Above  20  psi,  the  rates  appeared  to  continue  to 
increase  identically.  The  data  in  Figure  11  are  admittedly  too  limited 
to  warrant  a  positive  confirmation,  but  like  the  data  in  Figure  10, 
they  c.croainly  tend  to  corroborate  the  loading-rate  hypothesis. 

Under  field  loadings,  stresses  in  the  superseismic  region  probably 
rise  more  uniformly  to  peak  intensity,  i.e.,  there  is  no  slow  lead-in. 
Therefore,  because  of  the  relatively  slow  lead-in  characteristics  of  the 
laboratory  loading  pulses,  the  UX  test  specimens  may  not  have  been 
loaded  fast  enough  initially  to  be  representative  of  field  conditions, 
and  hence,  the  seismic-associated  stress  le’ els  assumed  for  the  pretest 
property  recommendations  would  have  been  too  low.  From  the  above 
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reexamination  of  the  available  laboratory  data,  it  now  appears  that  a 
lower  bound  value  of  50  psi  is  more  reasonable  for  the  stress  level  as¬ 
sociated  with  seismic  velocity.  When  the  properties  derived  from  the 
AFK  analysis  of  the  CIST  field  data  become  available,  they  may  give  a 
better  indication  of  the  actual  maximum  stress  level  associated  with  the 
seismic  speeds.  However,  it  is  possible  that  CIST-derived  data  might  be 
an  upper  limit,  since  the  loading  rates  in  that  experiment  are  consider¬ 
ably  faster  than  those  in  Event  III,  except  perhaps  in  the  immediate 
vicinity  of  GZ. 

U.3  EFFECT  OF  HORIZONTAL-TO-VERTICAL  ANISOTROPY 

It  was  obvious  from  the  presence  of  bedding  planes  that  response  ol 
the  Kayenta  rocks  would  be  governed  to  some  degree  bj-  anisotropic  prop¬ 
erties.  Data  from  direct  pull  and  Brazil  tension  tests  confirmed  that 
these  materials  have  a  substantially  higher  tensile  strength  in  the 
horizontal  direction  than  in  the  vertical  direction  (Reference  l).  Un¬ 
confined  compression  (UC)  tests  were  conducted  by  Terra  Tek,  Inc.  (TT), 
on  core  with  bedding  planes  oriented  parallel  and  perpendicular  to  the 
loading  axis  (Reference  10 ).  Strengths  from  the  parallel  oriented  tests 
were  as  much  as  50  percent  lower  than  those  obtained  from  perpendicu¬ 
larly  oriented  tests. 

Recent  compressibility  data  from  WES  UX  tests  on  Kayenta  sandstone 
specimens.  «ome  of  which  contained  thin  clay  seams  or  beddings  and  some 
oi  which  were  uniform,  are  shown  in  Figure  12.  As  expected,  the  uni¬ 
form  specimens  were  less  compressible;  their  responses  to  horizontal  or 
vertical  loadings  were  essentially  the  same.  On  the  other  hand,  verti¬ 
cally  oriented  specimens  containing  horizontal  clay  seams  were  much  more 
compressible  than  the  horizontally  oriented  test  specimens.  Since  the 
tests  were  limited  to  the  relatively  competent  unweathered  sandstones  of 
the  Kayenua  formation,  nothing  is  known  as  yet  about  the  anisotropic 
compressibility  characteristics  of  other  materials,  such  as  the  weath¬ 
ered  siltstone. 

As  previously  indicated,  the  codes  used  in  the  p.etest  ground 
motion  prediction  calculations  employed  isotropic  material  models . 
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Therefore,  even  though  the  laboratory  data  obtained  for  the  Mixed  Com¬ 
pany  sandstones  definitely  indicate  anisotropic  stress-strain  and  ten¬ 
sile  strength  characteristics,  a  basic  problem  still  remains  as  to  how 
to  utilize  this  information  in  specifying  "effective"  isotropic  property 
for  subsequent  calculations. 

4. It  REClNT  IN  SITU  AND  LABORATORY  STRENGTH  DATA 

One  of  the  pretest  analysis  assumptions  was  that  the  uhen-available 
laboratory  test  data  on  Kayenta  rock  specimens  were  representative  of  in 
situ  properties,  since  the  Rock  Quality  Index  of  the  core  was  greater 
than  90  percent  and  no  other  data  existed  that  could  be  used  to  adjust 
the  laboratory  data  to  reflect  in  situ  conditions.  More  recently,  how¬ 
ever,  a  UC  test  was  conducted  by  IT  at  a  rock  outcrop  near  the  CIST  ex¬ 
periment  on  a  large  block  of  in-place  sandstone  (Reference  10).  Corres¬ 
ponding  laboratory  tests  were  also  conducted  on  this  sandstone.  The 
field  strength  was  approximately  30  percent  lower  than  the  strength  from 
similarly  oriented,  small-size  laboratory  specimens  of  the  same  sand¬ 
stone.  In  the  absence  of  comparative  data,  on  confined  specimens,  it  is 
assumed  that  the  unconfined  results  represent  an  upper-bound  reduction 
factor,  i.e.,  in  situ  and  laboratory  test  data  should  agree  best  when 
the  materials  are  confined  under  large  static  overburden  or  dynamic  live 
loadings . 

Additional  information  supplied  by  Lawrence  Livermore  Laboratory 
(LLL)  on  one  type  of  Mixed  Company  sandstone  has  direct  bearing  on  the 
validity  of  the  pretest  failure  envelope  shapes  at  high  pressures  (Ref¬ 
erence  ll).  Prior  bo  the  LLL  tests,  the  only  failure  data  available 
were  the  WES  low-pressure  results  and  those  obtained  by  TT  at  confining 
pressures  of  2  and  4  kbars  (see  solid  circles  in  Figure  13)-  A  nearly 
linear  relationship  was  assumed  as  shown  in  the  figure  in  order  to  ex¬ 
tend  the  envelope  beyond  that  range  of  pressure.  The  LLL  data  confirmed 
the  earlier  WES  and  TT  results  but  not  the  assumed  high-pressiire  enve¬ 
lope.  Rather,  they  showed  an  extremely  nonlinear  envelope  that  sub¬ 
stantially  increased  under  confining  pressures  in  excess  of  4  kbars 
(see  open  circles  in  Figure  13). 
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LLL  also  determined  the  failure  strengths  of  specimens  that  had 
been  first  hydrostatically  loaded  to  7  kbars  and  then  unloaded  to  lower 
confining  pressures  prior  to  application  of  shearing  stresses.  Struc¬ 
tural  collapse  apparently  occurred  under  the  7-kbar  preloading  since 
the  resulting  failure  envelope  was  nearly  linear,  similar  to  that  ex¬ 
pected  for  a  dense  cand  (see  crosses  in  Figure  13).  Thus,  the  sand¬ 
stone  material  could  have  two  failure  envelopes,  depending  on  its  load¬ 
ing  history.  Although  such  crushing  could  only  occur  under  very  intense 
pressures  of  the  magnitude  expected  at  ranges  very  close  to  GZ,  it  could 
affect  the  characteristics  of  wave  forms  propagating  out  from  that 
region. 

4.5  RECOMMENDED  POSTSHOT  CONSTITUTIVE  PROPERTIES 

The  complete  set  of  representative  stress-strain  and  strength  rela¬ 
tions  for  Layers  17.  through  7  is  presented  in  Appendix  A.  They  consist 
of  UX  axial  stress-axial  strain  relations,  UX  stress  paths  for  both 
loading  and  unloading,  and  TX  failure  envelopes.  The  recommended  post¬ 
shot  profile  and  composition  properties  were  presented  in  Table  1,  and 
the  revised  seismic  velocity  values  were  given  in  Table  2. 

As  stated  in  Section  4.2,  initial  UX  compression  moduli  were  recal¬ 
culated  for  all  of  the  site  materials  using  the  revised  seismic  velocity 
values.  Each  of  the  UX  relations  for  Layers  1A  through  7  reflects  these 
moduli  (Appendix  A).  The  constitutive  properties  for  the  two  overburden 
layers  reflect  the  effect  of  increased  water  in  the  surface  layer; 

Layer  1A  is  initially  more  compressible,  locks  at  a  smaller  axial  strain, 
and  has  less  strength  than  Layer  IB. 

For  the  Kayenta  materials  (Layers  2,  3,  and  4),  the  initial  por¬ 
tions  of  the  originally  recommended  UX  curves  were  altered  to  reflect  a 
50-psi  stress  level  associated  with  their  seismic  velocities.  The  new 
curves  were  then  smoothly  transitioned  so  as  to  merge  into  the  previ¬ 
ously  recommended  relations  at  a  stress  level  of  about  2,000  psi,  such 
as  shown  in  Figure  l4.  No  changes  were  recommended  for  stress  levels 
above  2,000  psi,  primarily  because  of  the  uncertainty  expressed  in  Sec¬ 
tion  4.3  as  to  how  to  utilize  the  various  horizontal  and  vertical  data 
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to  specify  "effective"  or  "average"  isotropic  property  values. 

The  unconfined  strength  data  for  Layer  2  were  reduced  by  30  percent 
from  the  pretest  envelope  to  account  for  differences  between  in  situ  and 
laboratory  test  conditions.  The  envelope  was  gradually  merged  back  to¬ 
ward  the  values  of  the  pretest  envelope  at  higher  mean  normal  stress 
levels.  The  principal  stress  difference  levels  associated  with  uncon¬ 
fined  strength  for  Layers  3  and  were  reduced  by  20  percent  for  the 
same  reason,  and  the  yield  surfaces  were  then  merged  back  into  the 
values  of  the  pretest  envelopes.  As  shown  in  Figure  15,  such  changes 
are  almost  trivial  for  envelopes  depicting  strong  Coulomb-type  initial 
response.  More  substantial  changes  were  made  at  high  stress  levels 
(20  to  60  ksi  mean  normal  stress)  where  the  envelopes  reflect  the  non¬ 
linear  shape  previously  observed  by  LLL  (see  Figure  13). 

The  posttest  constitutive  properties  for  tee  Wingate  Formation, 

Layer  5s  are  the  same  as  the  pretest  properties  with  the  exception  that 
the  initial  portion  of  the  UX  curve  has  been  modified  to  reflect  the 
50-psi  seismic  precursor.  The  remaining  Chinle  and  Precambnan  mate¬ 
rials  (Layers  6  and  7)  are  assumed  to  be  elastic,  with  both  sets  of  prop¬ 
erties  reflecting  the  revised  compression  wave  velocities  of  11,000  and 
18,000  ft/sec,  respectively. 


Figure  8  Recommended  UX  relations  for  Layers  1A,  IB,  f’ud 
alternate  Layer  1C. 
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Figure  10  Comparison  of  loading  times  versus  peak  stress  for  laboratory  tests  with  those 
from  field  test  measurements . 
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Figure  12  Results  of  horizontally  and  vertically  oriented  UX  tests 
on  uniform  Kayenta  sandstone  specimens  and  specimens  with  clay  seams. 


kl 


k 

kS 

13 1 

kC3 
to  o 

kk 


\ 

X\ 


syvaM  £Jd-zd 

ddld  SS3H19  IVdIDNIHd 


Figure  13  Preliminary  results  of  TX  tests  by  LLL  on  virgin  specimens  of  Mixed  Company  sandstone 
(circles)  and  specimens  previously  subjected  to  7~hbar  hydrostatic  confinement  (crosses). 
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Figure  lH  Comparison  of  UX  loading  and  unloading  relation  recommended 
for  posttest  Layer  3  (solid  lines)  with  relation  recommended  for  pre¬ 
test  Layer  III  (dashed  lines). 
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Figure  15  Comparison  of  low-pressure  TX  failure  envelope  recommended 
for  posttest  Layer  3  (solid  line)  with  relation  recommended  for  pretest 
Layer  III  (dashed  line). 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  pretest  ground  motion  predictions  did  not  agree  very  well  with 
the  ground  shock  measurements  made  during  Mixed  Company  Event  III.  How¬ 
ever,  a  number  of  invalidated  assumptions  were  made  in  the  process  of 
developing  the  preshot  profile  and  property  idealizations.  ?.n  addition, 
the  mathematical  models  were  hurriedly  fit  in  order  to  make  the  pretest 
calculation  and  did  not  closely  replicate  all  features  of  the  recom¬ 
mended  profile  and  properties.  As  a  result,  it  was  concluded  that 
sufficient  uncertainties  existed  to  warrant:  (l)  a  reevaluation  of  the 
profile  and  properties  in  the  light  of  more  recent  data,  and  (2)  a  re¬ 
calculation  of  Event  III  with  a  much  more  precise  fit  to  the  constitu¬ 
tive  properties  and  profile  resulting  from  this  evaluation. 

The  pretest  calculational  grid  extended  only  to  a  depth  of  about 
300  feet,  i.e.,  into  the  Wingate  Formation.  Examination  of  possible 
travel  paths  showed  that  signals  reflected  from  the  Chinle  and  Precam- 
brian  layers  could  have  arrived  in  time  to  influence  the  ground  motions 
measured  during  Event  III.  It  is  therefore  strongly  recommended  that 
the  grid  used  in  any  recalculation  extend  to  depths  of  500  to  600  feet, 
i.e.,  into  the  Precambrian  basement. 

Samples  taken  just  a  few  days  after  the  event  showed  that  the  wet 
weather  conditions  existing  at  shot  time  resulted  in  a  significant  in¬ 
crease  in  the  water  content  of  the  upper  2  to  3  feet  of  the  overburden 
soil  over  that  determined  earlier;  tests  on  these  samples  showed  that 
the  increased  water  substantially  affected  the  compressibility  of  these 
materials.  It  is  now  recommended  that  the  postshot  calculational  zoning 
be  set  fine  enough  to  accommodate  an  idealization  of  the  overburden 
soil  into  two  horizontal  layers,  rather  than  into  one  as  orginally  rec¬ 
ommended  for  preshot  calculations. 

With  the  uncertainties  surrounding  the  definition  of  seismic  veloc¬ 
ities  and  constitutive  properties,  it  appears  that  the  most  practical 
approach  at  this  time  would  be  to  keep  the  profile  for  the  planned  re¬ 
calculation  as  simple  as  possible.  Thus,  it  is  recommended  that  the 


U-foot-thirk  artificial  layer,  which  was  included  in  the  preshot  profile 
to  represent  ;he  soft  clayey  conglomerate  material  occurring  randomly 
throughout  the  Kayenta  Formation,  he  eliminated  from  the  postshot  pro¬ 
file.  However,  the  basic  assumption  of  lumping  small  zones  of  material, 
which  cannot  be  incorporated  in  the  profile  due  to  computational  grid- 
size  limitations,  into  one  or  more  larger  layers  should  be  evaluated  in 
future  parametric  jcudies. 

An  extensive  reanalysis  of  the  available  surface  refraction  survey 
data  along  with  analysis  of  the  Event  III  ground  motion  data  resulted  in 
revisions  to  the  seismic  velocity  profile  as  listed  in  Table  2.  Hori¬ 
zontal  velocities  measured  by  AFWL  in  its  CIST  experiment  could  not  be 
reconciled  with  the  refraction  survey  data;  a  cross-hole  seismic  survey 
is  definitely  recommended  to  assist  in  resolving  this  uncertainty. 

Based  on  the  data  and  analyses  presented  herein,  new  constitutive 
property  recommendations  have  been  developed  for  each  of  the  layers  in 
the  proposed  postshot  profile  (Appendix  A).  Values  of  initial  UX  moduli 
were  recalculated  for  all  of  the  site  materials  using  the  revised  seis¬ 
mic  velocity  values.  The  stress  levels  associated  with  these  moduli 
were  originally  assumed  to  be  quite  low,  e.g.,  less  than  10  psi  for  the 
Kayenta  Formation  materials.  A  lower  bound  value  of  50  psi  is  now  rec¬ 
ommended  as  being  more  reasonable  for  the  initial  Kayenta  "precursor" 
stress.  This  phenomenon  appears  to  be  a  function  of  loading  rate  and 
definitely  requires  further  study. 

Tests  on  the  Kayenta  materials  also  revealed  that  horizontal-to- 
vertical  anisotropy  undoubtedly  affected  the  ground  motion  results.  But 
the  problem  of  how  to  utilize  the  various  horizontal  and  vertical  data 
to  specify  meaningful  "effective"  or  "average"  isotropic  property  values 
still  remains  as  a  nagging  item  for  x’urther  research. 

In  the  preshot  analysis,  yield  envelopes  were  extended  in  rela¬ 
tively  simple  linear  fashion  to  pressure  ranges  beyond  that  of  the  pre¬ 
test  data.  As  a  xesult  of  recent  LLL  data,  substantial  changes  have 
been  made  in  the  failure  envelopes,  which  now  reflect  a  highly  nonlinear 
behavior,  including  significant  strength  increases  at  high  pressures. 

Care  snould  be  taken  in  fitting  the  material  models  for  the 
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proposed  postshot  calculation  in  order  to  insure  that  they  reflect  the 
revised  properties  as  closely  ar  possible.  Recent  developments  related 

I 

to  cap-type  constitutive  models  have  significantly  improved  their  abil-  j 

ity  to  mirror  detailed  and  complex  material  property  specifications.  If 
after  doing  this  the  recalculated  ground  motions  still  do  not  reasonably 
match  the  field  data,  attention  should  be  directed  torard  defining  prop¬ 
erties  and  developing  simplified  models  for  loading- rate-dependent 
and/or  anisotropic  materials. 
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APPENDIX  A 


REPRESENTATIVE  STRESS-STRAIN  AND  STRENGTH  RELATIONS 

An  idealized  profile  and  set  of  constitutive  properties  were  devel¬ 
oped  for  use  in  the  postshot  calculation  of  the  Mixed  Company  Event  III 
test.  The  recommended  postshot  profile  and  composition  properties  are 
summarized  in  Table  A.l  which  lists  the  layer  number,  the  material  de¬ 
scription,  the  depth  range  for  each  layer,  the  values  of  wet  density, 
water  content,  and  volume  of  air,  the  seismic  velocity,  and  the  tension 
limit  (tension  cutoff)  of  axial  stress.  The  constitutive  properties 
consist  of  a  UX  axial  stress-axial  strain  relation,  a  principal  stress 
difference-mean  normal  stress  path  for  UX,  and  a  TX  failure  envelope. 

The  constitutive  properties  for  Layer  1A  are  shown  in  Figures  A.l  to 
A. 4 ;  Layer  IB  properties,  in  Figures  A. 5  to  A. 10;  Layer  2,  in  Fig¬ 
ures  A. 11  to  A. 15;  Layer  3,  in  Figures  A.l6  to  A. 19;  Layer  4,  in  A. 21 
to  A. 23;  Layer  5,  in  A. 24  to  A. 26;  and  Layers  6  and  7,  in  Figure  A. 27. 
All  of  the  properties  reflect  the  response  of  the  materials  to  live 
stress  loadings.  Tension  cutoff  values,  including  the  contribution  due 
to  overburden  stress,  are  indicated  for  layer  interface  locations,  i.e., 
top  and  bottom  of  each  layer. 
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'ABLE  A.l  SUMMARY  OF  POSTSHOT  PROFILE  AND  COMPOSITION  PROPERTIES  FOR  MIXED  COMPANY  ANALYSES 


Representative  az  versus  ez  relation  for  uniaxial  strain  with  unloading  curves 
100 ,  500,  and  1,150  psi  for  Layer  1A. 
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Figure  A. 3  Representative  az  versus  ez  relation  for  uniaxial  strain 
to  oz  =  6,000  psi  for  Layer  1A. 
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Figure  A.  5  Representative  cfz  versus  ez  relation  for  uniaxial  strain 
with  unloading  curves  from  az  =  100  and  500  psi  for  Layer  IB. 
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Figure  A. 6  Representative  oz  versus  ez  relation  for  uniaxial  strain 
with  unloading  curves  from  oz  =  500  and  1,950  psi  for  Layer  IB. 
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Figure  A. 7  Representative  oz  versus  ez  relation  for  uniaxial  strain 
to  a7  =  30,000  psi  for  Layer  IB. 
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Figure  A. 12  Representative  az  versus  ez  relation  for  uniaxial 
strain  with  unloading  curves  from  crz  =  5")0  and  1,000  psi  for  Layer  2. 
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Figure  A. 13  Representative  az  versus  ez  relation  for  uniaxial 
strain  with  unloading  curves  from  az  =  500,  1,000,  5 ,000,  and 
8,000  psi  for  Layer  2. 


Representative  (az  -  ar)  versus  p  stress  paths  for  uniaxial  strain  and 
versus  p  failure  envelope  for  triaxial  shear  to  p  =  1,200  psi  for  Layer 


Representative  (oz  -  ar)  versus  p  stress  path  for  uniaxial  strain  and 
versus  p  failure  envelope  for  triaxial  shear  to  p  =‘6,000  psi  for  Layer 


AXIAL  STRESS 


Figure  A. IT  Representative  az  versus  ez  relation  for  uniaxial 
strain  with  unloading  curves  from  az  =  500,  1,000,  2,000,  and 
8,000  psi  for  Layer  3* 
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Figure  A.18  Representative  (oz  -  ar)  versus  p  stress  paths  for  uniaxial  strain  and 
(cz  -  ar)  versus  p  failure  envelope  for  triaxial  shear  to  p  =  1,200  psi  for  Layer 
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Figure  A.  22 
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Representative  cz  versus  ez  relation  for  uniaxial  strain  with  unloading 
oz  =  200  and  500  psi  for  Layer  5. 


Figure  A.  25  Representative  az  versus  ez  relation  for  uniaxial 
strain  with  unloading  curves  from  a  =  2,000  and  8,000  psi  for 
Layer  5. 
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